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We study the head-on collision of two equal mass, nonrotating black holes. We 
consider a range of cases from holes surrounded by a common horizon to holes initially 
separated by about 20M, where M is the mass of each hole. We determine the 
waveforms and energies radiated for both the t = 2 and i = 4 waves resulting from 
the collision. In all cases studied the normal modes of the final black hole dominate 
the spectrum. We also estimate analytically the total gravitational radiation emitted, 
taking into account the tidal heating of horizons using the membrane paradigm, and 
other effects. For the first time we are able to compare analytic calculations, black 
hole perturbation theory, and strong field, nonlinear numerical calculations for this 
problem, and we find excellent agreement. 
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I. Introduction. The collision of two black holes is considered to be one of the most 
promising and important astrophysical sources of detectable gravitational radiation in our 
Universe Since LIGO |I| and VIRGO [[[]] are expected to begin taking data during this 
decade, it is important to perform accurate calculations detailing the shape and strength 
of the waveforms generated during such events. The information gained from the detected 
waveforms should allow one to reconstruct the astrophysical parameters of the system, and 
will provide the first direct and unambiguous evidence for the existence of black holes if the 
unique signature of the quasinormal modes of the hole is excited. 

In this paper we present results of numerical studies of time symmetric, axisymmetric 
head-on collisions of equal mass black holes. We have been able to extract the waveforms 
and the total energy emission resulting from the collision. Analysis of the waveforms reveals 
clearly for the first time that the quasinormal modes of the final black hole are strongly 
excited from the collision. This work extends and refines the early work of Smarr and 
Eppley |3|,|]] that suggested that the normal modes of the final hole were excited, but the 
resolution and waveform extraction techniques available at that time did not permit a clean 
and unambiguous matching to black hole normal modes. The numerical difficulties inherent 
in this problem also led to fairly large uncertainties in the total energy radiated (Smarr 
quotes a probable uncertainty factor of 2 Because of the importance of this fundamen- 
tal physical problem, we have revisited this calculation with the benefit of more powerful 
computers and improved analytic and numerical techniques developed over the intervening 
15 years to calculate unambiguous waveforms and energy fluxes resulting from the collision. 

The initial data sets that we adopt are analytic solutions originally discovered by Mis- 
ner || and subsequently analyzed and evolved by Smarr and Eppley. They are characterized 
by a parameter fi determining the mass M of a hole and the proper separation L/M of the 
two holes (see Table 1), and consist of two throats connecting identical asymptotically flat 
spacetimes In this paper we apply the code described in Ref. to compute evolutions for a 
family of Misner spacetimes representing equal mass black holes colliding from distances of 
between 4M and 20M. If the throats are close enough together (for /i < 1.36) [0] a common 
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apparent horizon surrounds them both, so that the system really represents a single, highly 
perturbed black hole. For throats separated by more than about 8M we are confident that 
there is no common event horizon surrounding them, as shown by directly integrating light 
rays (see also a hoop conjecture argument ||). 

II. Waveforms. Numerically generated waveform templates will be essential for analysis 
of data collected by future gravitational wave detectors U . The method we use to calculate 
waveforms is based on the gauge invariant extraction technique developed by Abrahams and 



Evans [0 and applied in Ref. JL1]] to black hole spacetimes. 

For all of the cases studied in this paper we have extracted both the i = 2 and £ = 4 
waveforms at radii of 30, 40, 50, 60, and 70M. By comparing results at each of these radii we 
are able to check the propagation of waves and the consistency of our energy calculations. 
In general we find agreement among detectors at different radii to be within 10-20%. We 
have also run all the calculations with resolutions of 200x35 (radial, angular) and 300x55 
with only minor differences in the 1 = 1 waveforms and the total energy radiated. (But as 
detailed in Ref. |J the amplitude and precursor of the I = 4 waveform are rather sensitive 
to computational parameters.) 

In Fig. 1 we show the I = 2 extracted waveform for the case // = 2.2, for which there is 
no initial common event horizon. The solid line shows the waveform detected at a distance 
r = 40M and the long dashed line shows the waveform extracted at r = 60M. The wave 
is clearly propagating away from the hole at light speed with essentially invariant shape 
and amplitude, with a wavelength of 2 x 16. 8M, confirming the original findings of Smarr 
and Eppley 0. However, our more accurate code now allows us to go beyond estimating 
the wavelength and to fit quantitatively the waveform to results known from black hole 
perturbation theory. The short dashed line shows the result of fitting the r = 40M waveform 
(in the range 64M < t < 160M) to a linear combination of the fundamental and first excited 
state of the I = 2 quasinormal mode for the final black hole with mass 2M. The fit is quite 
good, matching both the wavelength and damping time, showing that the final black hole 
mass is indeed very close to the total mass of the spacetime. 



In Fig. 2 we show the I = 4 waveform for the same case, extracted at the same radius. 
Again this waveform has been fit (over a similar range) to a superposition of the fundamental 
and first excited I = 4 quasinormal modes of the black hole. Waveforms for all cases studied 
show similar behavior: the normal modes of the final black hole are excited and account for 
most of the emitted signal. 

III. Energy Output. We first present a semi-analytic determination of the energy radiated 
by the colliding black holes. Such an estimate is interesting not just as a confirmation of the 
numerical results, but more importantly, it provides physical understanding of the numerical 
data. 

III. A. Analytic Estimate. We base our understanding of the two black hole collision on 
the well-studied problem |12|JT3| of a test point particle m at rest from infinity plunging into 



a Schwarzschild black hole M. The gravitational wave energy radiated for the test particle 



case is found to be 12 



E = 0.0104m 2 /M ,/orm« M. (1) 

We build on this accurate result to approximate the energy loss for the case of two unequal 
mass black holes of mass m and M, with m < M. We include correction factors into Eq. 
(H), so that it can better approximate the two black hole collision. The most important 
effects one has to take into account |J are that, (i) m may be comparable to M, (ii) the 
infall is not from infinity, (Hi) the black hole, unlike a point particle, has internal dynamics. 

It is useful to first understand the m 2 /M dependence of Eq. ([!]). In Newtonian approxi- 
mations, the quadrupole moment of the system is / ~ mr 2 for m <C M, where r is the radial 
distance between m and M. The gravitational wave luminosity L is given by 

Loci 2 ~ m 2 (rrf ~ m 2 (V 3 /r 5 ) . (2) 

As most of the radiated energy is emitted when m is falling near the horizon of M, to obtain 
the total energy radiated E, we can evaluate the luminosity at r ~ M and multiply it by a 
time scale of M. This gives E cx m 2 /M as in Eq. ([!]). Now for m not much smaller than 



M, it is more accurate to use / ~ fir 2 , with /i = (mM) / (M + m). The quadrupole moment 
formula then leads to ]|] E oc /i 2 /M, and we modify Eq. ([[]) to be 

E = OMOAff/M. (3) 

Notice that for m = M, the /i 2 of Eq. introduces a quite significant factor of 1/4. 

The quadrupole formula (2) also suggests how the expression should be modified when 
the infall is not from infinity. There are two effects: (i) there is less time to radiate when 
falling from a finite distance, and (ii) the infalling velocity is smaller for the same separation. 
The latter effect is much more important. These two effects can be combined into a single 
factor F ro , 

F ro = ErjEoo = / r(r) 2 dr/ / r(r) 2 dr. (4) 

Here r is the radial Schwarzschild coordinate, r Q is the starting point of the 
infall, and for the infall velocity f we use the relativistic expression (1 — 



2M/r)^2M/r - 2M/r /^Jl - 2M/r Q |Lj. We find F To ranges from ~ 0.4 (for r = 6M) to 
1 (for r Q = oo). 

All considerations up to this point are the same whether the infalling object is a point 
particle or a black hole. As far as the gravitational wave output is concerned, the most 
important difference between a point mass and a black hole is that a black hole has internal 
dynamics. There are more channels into which the initial gravitational potential energy 
in the system can dissipate. Such dissipation decreases the kinetic energy and hence the 
velocity of the infalling hole, reducing the wave output. 

When the hole m is falling in the static gravitational field generated by M, the horizon 
of m will be deformed by the tidal force. In the membrane paradigm ]K[ of a black hole, in 
which the horizon is treated as a 2-D surface living in a 3-D space, endowed with physical 
properties such as viscosity, this tidal deformation heats up the horizon. The heating is 
described by the horizon equations |T^JT^| — u a b + (g — G)a ab + {2a ac + ■jac^)^ = tab, —® + 



g6 — 9 2 /2 = <T a b<r ab , and 7^ = 2a a b + labO- Here jab is the 2-D metric of the horizon, 9 and 



a a h are the expansion rate and shear of the horizon generators, and g = l/(4m) is the surface 
gravity, while e a t, is the normalized electric part of the Weyl tensor [C a ^ u l^l v , with I 11 the 
horizon generators] , here to be evaluated on the horizon of the infalling hole m as it falls in 
the external tidal field [Mj r 3 ) of the hole M. For our present purpose, it suffices to solve 
these equations in the linear approximation, dropping all terms quadratic in the horizon 
deformation. The Green's function for the linearized equations satisfying the teleological 
boundary condition is simply G(t,t') = exp[g(t — t')] for t < t' and G(t,t') = for t > t'. 
Hence 9 can be obtained by a simple integration. The fraction of the available gravitational 
potential energy originally in the system which is dissipated into the heating of the horizon 
of the infalling hole is given by 

Am If 1 r 2m+2M 

The integration in Eq. [5] is terminated at the point when the two holes are engulfed by a 
common horizon. For m = M, the heating on the horizon of M is the same as that on m, 
so the total fraction of energy going into horizon heating is given by two times Eq. [5], i.e., 
the reduction factor for the energy available for wave generation is = 1 — 2Am/m, which 
decreases with increasing initial separation. For the range of r a studied here, Fh gives a 
reduction from about 3 to 13%. 

With these correction factors taken into account, the total gravitational energy radiated 
is given by 

E = 0.010V/ M x F vo x Fh- (6) 

There are other correction factors in Eq. (0) that we studied and included in Ref. |J. As 
their effects are much smaller, we shall not discuss them here. We note that for holes initially 
infinitely separated, F To is 1, while Fh is 0.86. In the next subsection we compare Eq. (^) 
for the total energy radiated to numerical results. 

III. B. Numerical Results. The waveforms shown in Figs. 1,2 have been normalized so 
that the total energy carried away from the system is given by E = ip 2 / (32n) for all I modes, 
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where ip is the extracted waveform. The total energy E is computed from this expression. For 
a check on the accuracy and consistency of this calculation, we have computed E at all five 
detectors listed above, and at both high and medium grid resolution. We have also computed 
the energy loss via two independent curvature based indicators, the Newman-Penrose Weyl 
tensor component ^4 and the Bel-Robinson poynting vector These constructions yield results 
consistent with the gauge invariant formulations when computed in the asymptotic far field 
spacetime |J. 

In Fig. 3, we compare results of various analytic estimates of the energy loss to our 
numerical results for the dominant i = 2 waveform calculation. The connected circles show 
the maximum possible radiation output obtained by comparing the initial black hole masses 
estimated by the areas of the horizons (or a single horizon if the holes are close enough) to the 
total mass of the spacetime. For large separation this number approaches 29% as expected 



from the work of Hawking |17 |. The six clusters of unconnected symbols result from our 
numerical simulations for the six cases studied (listed in Table I). Each symbol corresponds 
to the £ = 2 energy computed at each of five "wave detectors" located throughout the 
radiation zone as before. 

The analytic estimate for total gravitational wave energy output given by Eq. |6] is plotted 
as a dashed line. The agreement of the analytic and the numerical result is highly remarkable, 
as the analytic and numerical results were obtained in a double blind manner. For L/M 
less than about 9 the analytic formula overestimates the actual energy output computed 
numerically. This is to be expected since for small enough separation the holes are initially 
engulfed by a common event horizon and the approximation for colliding black holes is 
inappropriate. For reference, the early results of Smarr and Eppley are plotted as large 
crosses with error bars suggested by Smarr [|J. Within the large errors quoted, those early 
results are remarkably consistent with our results. 

III. Conclusions. We have performed numerical and analytic calculations predicting the 
gravitational waveforms generated and total gravitational wave energy emitted when two 
equal mass black holes collide head on. Both the I = 2 and i = 4 waveforms are fit 



remarkably well by a superposition of the fundamental and first excited state of the black 
hole quasinormal modes, and the total energy output is found to be on the order of 0.002M. 
The analytic study confirms and elucidates the numerical results, matching the numerical 
results remarkably well. The total energy output is far below estimates given by a simple 
application of the area theorem. Taken together, the analytic and numerical results indicate 
that even for holes that are initially infinitely separated, the total energy output will be the 
same order of magnitude. 
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for integrating light rays through the spacetimes, and Mark Bajuk for his work on visual- 
izations of our numerical simulations that aided greatly in their interpretation. This work 
was supported by NCSA, NSF Grant 91-16682, and NSERC Grant No. OGP-121857, and 
calculations were performed at NCSA and the Pittsburgh Supercomputing Center. 
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FIGURES 

FIG. 1. t = 2 waveforms for the case \i = 2.2. The solid line shows the waveform extracted at 
R = 40M and the long dashed line shows the waveform at R = 60M. The short dashed line shows 
the quasinormal mode fit. 

FIG. 2. £ = 4 waveforms for the case fi = 2.2. The solid line shows the waveform extracted at 
R = 40M and the short dashed line shows the quasinormal mode fit. 

FIG. 3. The total gravitational wave energy output E shown for various cases. The connected 
circles are the upper limit based on the area theorem, the clustered symbols show numerical results 
at various detector locations, the dashed line is the semi-analytic estimate, and the crosses show 
early results by Smarr and Eppley with their approximate error bars. 
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TABLES 





M 


L/M 


Apparent horizon 


1.2 


1.85 


4.46 


global 


1.8 


0.81 


6.76 


separate 


2.2 


0.50 


8.92 


separate 


2.7 


0.29 


12.7 


separate 


3.0 


0.21 


15.8 


separate 


3.25 


0.16 


19.1 


separate 



TABLE I. The physical parameters of the six initial data sets studied are summarized. M is 
the mass parameter defined in the text, L/M is the proper distance between the throats, and we 
note whether or not one apparent horizon surrounds both holes. 
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